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and we derive general amplitudes of spatially and temporally long ranged 
fluctuations in the considered diffusive system. 
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1. Introduction 

The description of dynamical processes in condensed matter greatly simplifies if 
fluctuations are studied which are slow and smooth compared to the microscopic length 
and time scales of the system. Then hydrodynamic equations for a small number of 
fields can be derived, either using rather general phenomenological considerations, or 
by coarse graining starting from a microscopic statistical mechanics description. In 
the latter a large number TV of particles needs to be handled and the hydrodynamic 
fields (normally) arise as coarse grained densities of conserved variables [Q. The 
Zwanzig-Mori operator formalism enables one to perform the coarse graining of the 
microscopic equations of motion using spatially Fourier transformed variables in the 
limit of vanishing wavevector q (corresponding to large wavelengths 2t: /q). Already in 
1931, Onsager explained how the microscopic equations in the limit oiq-^ Q and small 
frequency lead to the hydrodynamic equations j^]. He suggested that in an infinite 
system a perturbation described by macroscopic hydrodynamic equations decays from 
its inital value according to the identical dynamical equations as a long-wavelength 
and small frequency fluctuation around local thermodynamic equilibrium As an 
important side-product, this correlation functions approach has led to general and 
exact microscopic expressions for the phenomenological transport coefficients of the 
hydrodynamic equations (the Green-Kubo relations). 

Hydrodynamic equations, which are partial differential equations, require 
temporal and spatial boundary conditions to give unique solutions; see e.g. the 
examples in Q. Following Onsager, only the former are understood from microscopic 
many-body approaches, while, by studying infinite systems, the latter have been 
neglected. Within the phenomenological approach, simple continuity considerations 
lead to the required conditions on surfaces, yet their derivation from information 
about microscopic interactions and molecular parameters still appears desirable. 
First, this would provide rigorous statistical mechanics definitions of the parameters 
characterising the boundary condition; second, different conditions (like stick or slip 
for fluid flow, or flux versus no-flux with or without adsorption of particles at a 
surface) could be predicted from molecular interactions; and third, generalizations 
beyond the true hydrodynamic limit (e.g. for finite geometries) would become possible. 
Only the question about the tangential velocity of a flow along a solid boundary has 
a long history, which goes back to Maxwell, and, for rarefied gases, is quite well 
answered in the framework of the Boltzmann equation (Knudsen-layer problem ||5[]). 
Yet, beyond the dilute limit a fluctuating hydrodynamics calculation by Wolynes 
has uncovered the subleties arising from back flow patterns (coupling of hydrodynamic 
modes), and only rather recently there has been the flrst study of this problem in the 
microscopic correlation functions approach by Bocquet and Barrat 0|. In our study 
of the simpler system of a single conserved variable which macroscopically obeys a 
diffusion equation, we follow the approach of Bocquet and Barrat and connect a 
microscopic linear response calculation to the macroscopic hydrodynamic description 
via a generalization of Onsager's regression hypothesis. 

On the macroscopic level, the number density n( r, t) of interacting Brownian 
particles at the space point r and at time t obeys a diffusion equation 

dtn{v,t)^DV^ n{r,t) , (1) 

where D is the (gradient) diffusion coefficient which enters Eq. (|l|) as a 
phenomenological transport coefficient; dt = d/dt denotes a partial time derivative 
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and V — d/dr a spatial gradient. If the diffusing particles border a solid surface 
which moves with velocity v(t) then the number of particles (per unit area and time) 
displaced by the surface, n v, needs to be balanced by a particle flux j away from the 
boundary. As the latter obeys j = —DVn, the no-influx boundary condition on the 
solid surface becomes Q : 

e„ [n{ r, t) v{t) + D Vn( r, = , (2) 

here e„ is a unit vector normal to the surface, whose position is abbreviated as "bd" 
for boundary. 

In sections 2. and 3. of this manuscript, Eq. (||) will be derived up to linear 
order in v by coarse graining the many-body statistical mechanics description of 
interacting Brownian particles. At first the appropriate microscopic kernel is found 
(sect. 2.), and then its small q and oj behavior discussed (sect. 3.). The calculation 
entails the (standard) derivation of Eq. (^ including the Green-Kubo type calculation 
of D. The solution of Eqs. (|l|,§) around a spherical object to linear order in the 
perturbing velocity v is summarized in Appendix A for comparison reasons, while 
Appendix B contains technical material. Section 4. describes an application of our 
results. The power-law decay of the force experienced by a large sphere moving 
among the Brownian particles is deduced along with its mean-squared displacement, 
which exhibits a long time tail. 



2. Microscopic approach 

2.1. Smoluchowski equation, and notation 

The statistical mechanics basis for interacting Brownian particles is given by the 
Smoluchowski equation, which is a generalized diffusion equation in high-dimensional 
phase space |8[ ^ . It describes the temporal evolution of the many-body probability 
distribution ^'({ri},i), which depends on the positions r^ of all particles, i = 
1, . . . , N + 1, where we consider N identical bath particles with Brownian diffusion 
coefficient Di = Dq and one additional tracer with index .s = N + 1 and diffusion 
coefficient Ds'. 

^t^ = Y,D^^,■{^^- Fi) . (3) 

i 

Here di = d/drt, and energies are measured in units of the thermal energy. The 
particle interactions enter Eq. (^) via the potential forces = —djV{{ r^}) resulting 
from interactions between the bath particles (V^) and between tracer and bath 
particles (V), V — + V. Dynamic effects due to the background medium 
(hydrodynamic interactions) are neglected at the present stage. To reach the 
hydrodynamic limit we will take the size of the tracer to become much larger than 
the bath particle size. For the presentation in the main text, the tracer is assumed to 
be immobile from the outset, Ds = 0, and Appendix B verifies that the limit Z?s — > 
can be taken after the formal manipulations. Summation and indices are from now 
on always understood to run from 1 to iV, i.e., to exclude the index s = + 1 for 
the tracer. To simplify the presentation, we also introduce the radii as and a of the 
tracer and the bath particles, respectively. It is however important to realize that we 
consider arbitrary isotropic short ranged particle interactions, where a and may 
be effective state dependent sizes as e.g. in the case of soft repulsions of the form 
yp(r = |r, - rj-|)cxr-i2. 
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For the following, we introduce some further notational conventions. It is 
convenient to work with the backward or adjoint Smoluchowski operator 

n = DoY.{d, + F,) ■ , (4) 

i 

which gives the time evolution of variables A({ri}) on phase space: dtA — il.A. It 
also determines the time evolution of correlation (fluctuation) fmictions ^AB{t) ~ 
{A* exp {il.t} B) / {A* B) , which we normalize by their equal time values calculated by 
canonical averaging with the Gibbs-Boltzmann weight (...) cx / Jlila^ ■ ■ ■ ■ 
Note that here the tracer-particle interactions enter in full non-linear detail, and that 
the equilibrium weight is a stationary solution of Eq. (^). 

The fluctuating microscopic bath particle density at position r is given by 
q{v) — ^ - (5(r — Tj), with spatial Fourier transform, = ^ - expjiqrj}, where 
the g = contribution from the constant bulk density n will be neglected. The 
corresponding tracer density fluctuation is given by (f^ — expjiqrs}. Temporal 
Fourier decomposition shall be denoted by A^ = dt e*"*A(t), while the Laplace 
transformation is used with the convention: yl(a;) = dt e"^*A(f). 

2.2. Generalized Onsager regression hypothesis 

The connection between the statistical mechanics description on the Smoluchowski 
level and the macroscopic hydrodynamic picture shall be made by comparing the 
density fluctuations predicted from both descriptions for the identical given boundary 
problem in a simple geometry. 

In order to use the familiar Smoluchowski operator of Eq. (^, we consider the 
motion of particles around a spherical object, the tracer. Bath particle j experiences 
the short-ranged force Fj* — —djV close to it. Moving the tracer, by unspecifled 
external means, with vefocity v(t) induces a particle flux at its surface, which in 
linear order in v equals j'^'^(i) — n\-{t) on the macroscopic level. Deviations of the 
coarse grained particle density around the tracer n{r,t) from the bulk value n would 
enter in higher order in v only. The disturbance on the microscopic level, required 
to induce this applied particle flux thus can be obtained from requiring the non- 
equilibrium average of the tracer velocity to agree with the macroscopic value up to 
non-linear corrections: 

{dt r,)(-^ = v(0+O(v2) . (5) 

Adiabatically turning on the applied velocity in the infinite past eliminates initial 
value contributions in the deviatoric density, dn{r,t) = n(r, t) — n, and allows us to 
use a Fourier decomposition, v(i) — J (dijj/27r) e~''^* v^^. Linear response theory then 
connects the density deviation to the given disturbance via a (vector) susceptibility 
X(r, t). It vanishes for i < because of causality, and its spatial argument r is 
measured from the tracer sphere center. After Fourier transformation, both the 
macroscopic hydrodynamic result (cf. Appendix A) and the microscopic result (cf. 
sections 2.3 and 2.4) can be written as 

Sn^^^ ^ n • Xq(w) -HO(v^) . (6) 

Now, in the spirit of Onsager's hypothesis we assume that the microscopic 
calculation of Eq. (|6|) reduces to the macroscopic solution for smooth and slow 
fluctuations, viz. in the limit of small frequencies and wavevectors. Yet, in order 
to derive hydrodynamic boundary conditions, the coarse graining must be taken 
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with respect to the bath particle size a only, while the tracer size is required to 
satisfy Og S> a. Thus we keep fixed so that the macroscopic difi^usion equation 
description, while it applies for distances r ^ a only, nevertheless includes both far 
field (r ^ fls) and near field ( r <C Os). The latter case is equivalent to considering 
the density profile Sn{z,t) at a distance z close to a planar wall obtained formally 
when taking — > oo. Although this limit does not provide a faithful representation 
of a macroscopic boundary as an assembly of atoms, it has the virtue of being the 
conceptually simplest realization of a hydrodynamic boundary problem. 

Two aspects of the described approach are worth mentioning: First, while 
Eqs. (|^J^) are linear in the applied boundary flux, the particle-wall (tracer) 
interactions are included exactly. Thus on a local length scale the unperturbed 
equilibrium density variations arise, which somewhat differs from the approach to shear 
flow past a surface in Ref. . Second, as discussed e.g. by Kadanoff and Martin for the 
initial value problem ^ , a general perturbation to the fluid induces fluctuations in the 
non-conserved variables, which have to die out before the hydrodynamic description 
applies. For the present boundary perturbation the same reasoning applies, and thus 
the hydrodynamic description only holds for large distances, while locally deviations 
from the hydrodynamic solution need to appear; for rarefied gases these Knudsen- 
layers effects are familiar |^. In the present many-body linear response calculation 
the technical difficulty is connected to coarse graining across the equivalent layer, 
which has a width connected to the particle size a. 



2.3. Linear response calculation 

In order to proceed, the perturbation to the Smoluchowski operator needs to be 
found which gives the required velocity of the tracer in Eq. . Without hydrodynamic 
interactions, the perturbation equivalent to a constant solvent velocity, which is felt 
solely by the tracer, is by inspection 

An = v(t) • ds . (7) 

A standard linear response calculation using its adjoint Aft''' — — v(i) • ds, which acts 
on the probability density in Eq. (|^) ^J, gives the resulting deviation in an arbitrary 
variable A 

= = - f dT v(r)-(F, e"(*-") A) +0(v2) .(8) 

Thus, Eq. (^) becomes 

(9^ r,)("^) = ((r! + Ar!) r,)('-) = ( v(0 • 9. r,> + O(l?,/i?o, v') .(9) 

See Appendix B, for a more careful discussion for finite tracer diffusivities Ds > 0. 
As required, the perturbation Eq. (0) gives the average velocity of the tracer, which 
enters the macroscopic boundary condition Eq. (H). 

The linear response formula can also be applied to the microscopic density field 
g(r') at a (vector) distance r from the tracer center: r' — + r. Its unperturbed 
equilibrium value is proportional to the familiar tracer-particle pair correlation 
function ||l|: 5*(r) = {'^/n)^^{5[Y — {vi — r^)]), which gives the probability of 
finding bath particles at a distance r from the center of the tracer. The linear 
deviation in the density around the tracer induced by the perturbation, Eq. (R), follows 
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immediately from Eq. (^). and by comparison with Eq. the required hnear response 
susceptibility is found: 

n \l -\- luj ^ 

X( r, i) = - i ( F, e"* g{ r + r,)) e{t) , (10) 
n 

where the step- function 9{t) expresses causality. 
2.4- Time scale separation 

The exact linear response susceptibility varies on rapid microscopic time and length 
scales but also on smooth and slow ones, which are amenable to a hydrodynamic 
description. The Zwanzig-Mori projection operator formalism enables one to 
disentangle these contributions by splitting the resolvent into fast and slow subspaces 
[01 . The resolvent R{uj) = arises in the Laplace transform of a general correlation 
function, ^ab{oj) ~ {A* R{uj) B) / {A* B) , and contains poles which shift to vanishing 
frequency for smooth fluctuations {q ^ 0). These so-called hydrodynamic poles are 
connected with the exact conservation laws of the system, and the Zwanzig-Mori 
formalism isolates them. In the present situation, where Eq. (^) holds, there are 
only poles connected with particle number conservation: each particle, including the 
tracer, is conserved as is the total density, dtQq oc q for q ^ 0. While the more 
careful calculation in Appendix B takes into account the tracer, here we chose for 
the slow subspace the one spanned by the total density only. That is, we use the 
projector P = gq){NSq)~^{g*^, which is normalized by the equilibrium Brownian 
particle structure factor Sq = {g*qgq)/N The justification for this simplification 
is provided by the thermodynamic limit, in which only a non-extensive number of 
particles actually interacts with the tracer; see below and Appendix B. 

The exact identity obtained in the Zwanzig-Mori projection operator formalism 
| pO| gives for a general fluctuation function: 

{A*R{lj)B) = (A*R'{u;)B) + {A*{l+R'{Lj)n)PR{uj)P{l+nR'{uj))B) ,(11) 

where the reduced resolvent describes the fast dynamics decoupled from the slow 
fluctuations of the conserved density: 

R'i^)^Q nnn\- ^ with Q = 1 - P . (12) 

Thus, the coupling of the arbitrary variables A, B to the slow conserved density 
is found; explicitly it is obtained when writing out PR{uj)P in Eq. ([Tl|), and the 
slow variable couples in with static (i.e. {A*gc^) and frequency dependent (i.e. 
{A* R' {uo)Q,gc^) overlaps. 

In order to apply Eq. ( pi] ) to the correlation function in Eq. (p^, the (expected) 
problem arises that it is formed with variables that are not defined in a translationally 
invariant manner. Translational symmetry is broken by the boundary (i.e. measuring 
distances from the tracer). On the macroscopic level this could be handled by 
introducing the appropriate eigenfunctions that satisfy the boundary conditions for 
the prescribed geometry. Yet, on the microscopic level this would require determining 
the many-body eigenfunctions of the Smoluchowski operator Eq. (^) for a given 
force field arising from . Within the framework of fluctuating hydrodynamics 
Wolynes achieved a related task in a scattering-formalism calculation for the flow 
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of a Newtonian fluid past a wall Q. His calculation focused on the non-linear 
coupling of the hydrodynamic modes and thus could circumvent the study of the local 
variables close to the boundary. Consequently, he did not determine the boundary 
position microscopically and instead introduced a short distance cut-off (in his case 
irrelevant). Because we aim for an exact determination of the boundary condition, we 
chose the plane-wave decomposition of the density fluctuations which enables one to 
use Eq. (|l^), and apply it to the resolvent in a shifted coordinate system 

i?q(i) = g'^ i?(t) g% = R{t) (1 + OiDjDo)) . (13) 

It agrees with the original resolvent only if thermal tracer fluctuations are neglected 
(cf. Appendix B). In this limit, the Fourier transformed susceptibility becomes 

-n xq(a') = (F- R{u;) g^) = ( [l + i?'(^)r!] g^) <^> ^{lj) , (14) 

where we have introduced the tagged force density fluctuation F^ = F^fJq and the 
(normahzed) density correlator $q(a;) = {g*qR{uj)g q) /{NSq). Application of Eq. ( pl ) 
to the latter gives the familiar expression 



{g\ngq) (g*^nR'{L,)ngq) 



NSq NSg 



-l -1 


. ^ 2D0' 




-lUJ + q — 




Do _ 



The second expression in Eq. (|T5|) is taken in the hydrodynamic limit of small 
frequencies and wavevectors, where it gives the (transformed) fundamental solution of 
the diffusion equation Eq. (Q) . This leads to the known microscopic deflnition of the 
gradient diffusion coefflcient, D — Dq/Sq. Here, So is a normalized compressibility. 
The result for D may be called of Green-Kubo type because its apparent static nature 
originates in an instantaneously decaying associated current. 

The frequency independent (or instantaneous) overlap in Eq. (14) can be 
expressed in terms of the Fourier transform of the non-trivial part, ft,'*( r) — (7*( r) — 1, 
of the tracer-particle pair correlation function introduced above, (F'^pq) — zqrt/iq, 
but little further simpliflcation is possible in the retarded second term. Upon 
introducing the total force density fluctuation Fq = F^ e^'^'^^ , Q^g^ = 'iQ^j ^.^ 
Fq, and the final (still exact) result for x becomes (for t > Q) 



Xq{t)^-tqh%^qit)~t^ fdT{F^*R'{t-r) q- F< 

^ Jo 



W -(16) 



It is written as function of time to clearly present the instantaneous (first term) 
and retarded coupling of the density fluctuations to the susceptibility. Because of 
Newton's third law, the potential force Fs felt by the tracer can be reexpressed as 
the negative of the total force exerted by the tracer on all particles; F^ = —dsV = 
^ - cJi^* = — Fq. For the same reason it is opposite equal to the integrated total 
force; F^ = F^^g = — Fq^o = X^i^i^i = — Fq. A noteworthy aspect of the 
(straightforward) calculation in this section concerns the thermodynamic limit which 
is required in order for the obtained bulk quantities to take their standard values for 
an unbounded system. For example, the tracer bath particle interactions enter the 
expression for 5"^ (and consequently for D) via the equilibrium distribution function. 
Nevertheless, in the thermodynamic limit this correction vanishes because the assumed 
short-ranged interaction of the bath particles with the tracer decays beyond the 
distance of a few a, and the bulk of the particles is not affected. 
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3. Coarse graining and discussion 

The exact correlation function, Eq. (|l6|), describes the response of the system to 
an injected boundary flux of particles on all length scales. In order to derive the 
hydrodynamic boundary condition, coarse graining is required as discussed in section 
2.2. Appendix A collects the results from the macroscopic approach in order to 
compare them with the small wavevector and frequency limits of the microscopic 
susceptibility. 

3.1. Instantaneous response 

In order to familiarize oneself with Eq. (|l^), it is useful to consider a rapid velocity 
pulse on the tracer at time to, v(<) = V(5(t — t^), and to concentrate on the 
instantaneous response of the density: 

6n{r,t^ to) = nY ■ x{r,t = 0) = n V - V/i"(r) . 

This arises from the first term in Eq. ([l^), which simplifies because of $q(t = 0) = 1, 
and is determined by the equilibrium density profile of bath particles around the tracer. 
The inserted flux, nV, is packed close to the boundary according to the equilibrium 
fluid stucture /i* of the bath particles. The (Ursell) function h'^ varies between the 
universal limits, /i^ = — 1 for short distances where the hard core volumes of the 
particles are excluded by the tracer, and h'^ — far away from the tracer. In between, 
it shows layering over a distance of the order of a few a. 

In the hydrodynamic limit, which corresponds to Og ^ a here, the present work 
provides the connection of the position of the boundary to the molecular interaction 
potential V^. At a radial distance a from the tracer center, varies rapidly ^ and 
in the limit a —^ with fixed Og and r, may loosely be taken as a step function, 
/i'*(r) = —6{a — r). The macroscopic sphere asymptotically becomes impenetrable for 
the bath particles irrespective of the exact interaction potential. The latter however 
determines the exact boundary position cr, and its definition becomes: 

/iq -27rCT^ f{qa) + C(ct^) for > a and ga < I , (17) 

where /(x) = (sin a; — a;cosx)/a;^. Whenever Eq. ( p^) does not hold, possibly for long 
ranged forces or wetting situations, we expect Eq. (g) to be violated. Such situations 
are excluded in the following. A finite (positive or negative) surface excess density 
enters in the corrections of order cr^. In the following sections we show that exactly 
the same structure also appears in the retarded contributions to Eq. (p^, and that 
the boundary position a thus is a static equilibrium concept (see however Ref. for 
Newtonian fluid flow). 

In the limit of 3> a, the rapid variation of h'^ can be used to define a 
one-dimensional cut through the density profile, which in the limit ag ^ oo (and 
consequently ct ^ oo) would correspond to the situation at a wall jll. With the wall 
at a; = (T, and its normal vector pointing along x, the wall profile /i as function of 
r=r — (TX, a;=r-x follows 

h''{ v)~^h''^{x) + 0{f/a) for a ^ oo . (IS) 

It obeys, /i'*^(x — oo) ^ — 1 and h"^ {x oo) 0, with rapid variations on a 
length scale of order a around x ~ Q. Its one-dimensional Fourier transform is given 

by 

dx e'i-^ h'"^ {x) = ^ = —+H'{l + 0{q^a)),{l(i) 



Hydrodynamic boundary conditions 



9 



where the constant Ti.' is the surface density excess divided by n and is of order a itself. 
By shifting the origin to a; = {x = a), we ehminated (for simpUcity) a phase factor 
g»9i<^ in /ig^, which would prove convenient when keeping track of the wall position. 



3.2. Near-field solution 

Generically, boundary conditions are formulated when considering the motion in a half 
space bounded by a planar surface (wall). As discussed above for the instantaneous 
response, this situation can be realized in Eq. ( [l^ ) by taking the limit Og —> co and 
calculating x( r = crx + r, i) = i^^'^) ^ + /'^) to non-vanishing order. The 
result x^i^jt) describes the motion close to an infinite plane wall or, in general, the 
near-field solution for non-planar solid surfaces. Only its small wavevector limit is 
required in the following, and this simplifies because the force exerted on the diffusing 
particles by the wall (for — > oo) is perpendicular to it: 



1 + — (x- Fo R{uj) X- Fo 
n 



<f,.H- (20) 



Here we have used that for vanishing wavevector the reduced resolvent in the relaxation 
kernel again agrees with the full dynamics [||. The retardation kernel in Eq. ( pO| ) 
therefore has the familiar Green-Kubo form. If it could be replaced by a constant 
rate for small frequency, (x • Fq R{lu 0) x • Fq) T, then for consistency the 
square bracket would become [. . .] I + 0{q). Fortunately, in an exact calculation 
for vanishing concentration of hard Brownian spheres, n —^ 0, Dieterich and Pcschel 
have shown that 



( X . Fo Riuj) X • Fo) = J (1 + Oina^)) . (21) 
In this limit D = Dq, and the result of Eqs. (^,^ 20) agrees with the solution of the 



hydrodynamic equations, Eqs. (|l|,0), in the considered geometry; see Eq. ( A.3[ ). This 
proves the boundary condition Eq. (||) in the dilute limit of hard spheres 1 12[. 

For the general situation of interacting Brownian particles at finite concentrations, 
no exact calculations of the relaxation kernel in Eq. (|2C| ) are known. We proceed 
by performing a mode coupling approximation [p|,p^,[l4[, where the starting point 
is the more general expression of Eq. (|l^) as it captures near- and far-field 
terms. The conserved density fields are the slow variables and in the lowest pair- 
fluctuation approximation the overlap of the fluctuating forces with PkPk' needs to be 
considered.|^ In the small wavevector limit of interest, the memory function becomes 
identical to the well studied tracer force autocorrelation kernel, and its mode coupling 
result can be taken from the literature jl^ . 

^ (Fo R'it) q.Fo) « i^o q-hm (kk) ^^^^^ ^^(t) ,(22) 

where the tracer density fluctuation function does not appear explicitly (it equals 1 
because of Dg = 0). 

The mode coupling approximation of the relaxation kernel can be applied to 
the wall or near-field problem upon the realization that the forces arise from density 
fluctuations whose probability depends on the wall distance according to h'^^ (x) and is 

1 It can be expected that this approximation does not give numerically exact hydrodynamic 
results |l^ , but Schofield and Oppenheim jl^ have argued that this problem can be overcome by 
systematically taking the overlap with higher order density products into account. 
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independent of the paraUel coordinates, y and z. For the relevant wavevector region|^. 



the hi in Eq. (^22 



thus corresponds to 

(2^f <5,l(q) 



(23) 



where S\\ ( q) restricts the parallel wavevector to vanish, % = iz = 0, and the wall 
profile function H.'^ defined in Eq. (|l^) enters. We find the approximation: 

^ (i-Fo R{^) x.Fo) « {2nf 5||(q) '^V' ^° >(24) 



,27r 



where in the integrand ky 



0. Importantly, in the hydrodynamic limit w — > 0, 



the integral converges already for such small ^ 1/a that the structure functions 



can be replaced by their homogeneous zero wavevector limits, — > 1 and Sk 



So 



The latter is the bulk compressibility required to turn the single particle Dq into the 
gradient diffusion constant D. It is thus interesting that the exact one-dimensional 
result of Eq. (|2l] ) in the dilute limit only applies to the present case because the 
(isolated) particles experience no interactions and thus D = Dq in Eq. (|2|). Because 
the integration in Eq. (24) is dominated by fc ^ a, the density correlator can be 
replaced by its universal hydrodynamic limit from Eq. (^|). 

Collecting all terms together and performing the integrations gives the mode 
coupling approximation 



XdiLo)- x(27r)2 5||(q) <^ 1 + 



,(25) 



where the term in the square bracket arises from the considered pair density 
projections and . . . indicate higher order density projections which should be taken 
into account following an expansion procedure developed by Schofield and Oppenheim 
[p^ . Assuming the series of density projections in the square bracket in Eq. ( p5| ) to 
sum up to one, the comparison with Eq. (A. 3) proves the correctness of the boundary 
condition Eq. (0), now at finite concentrations. 



3.3. Far-field solution 

While the verification of the boundary condition Eq. (||) is achieved by the calculation 
of the near-field, it is instructive to also consider the density fluctuations very far 
away from the spherical boundary. In this limit, both particle sizes a and Og are small 
compared to the wavelength and the susceptibility in Eq. (f^ simplifies to 



Xq(0 ~ " with 



,s Do 
6n 



dt (Fo(O-Fo) ,(26) 



for qos ^ 1; where the density correlator takes its hydrodynamic form, ^q{t) = 
exp {^q^Dt) from Eq. (|l5|). The density profile around a moving tracer rearranges 
by particle diffusion and thus requires more and more time the larger the involved 
distances. In the steady case, w = 0, a power law density profile develops 



=o(r)^- 



(27) 

First, Eq. {||) is transformed to r-space and then h''{r) is analyzed according to Eq. jisl). The 
procedure is checked a posteriori from the convergence of the integrah 
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as follows from Eqs. (|6|,|15|,^6|) and transformation to r-space. It is of similar nature 
to the Oseen velocity profile around a colloidal particle in a Newtonian solvent Q , as 
such long-ranged patterns generally arise in hydrodynamic steady states [|l^. In the 
Brownian particle context, it is well known from calculations in the dilute limit, n — > 
[^,|l^,^, and the amphtude factor a extends those calculations to finite densities. 

Interestingly, the expression for a in Eq. ( p6| ) holds for arbitrary size ratios as /a, 
even beyond the macroscopic hydrodynamic limit, which is obtained for Os ^ a. For 
dilute hard spheres, the known result a{n — > 0, a/a^ — > 0) = — 27rcr^ (whe re a 



[P|,|l8| agrees with the expectation from the macroscopic calculation, Eq. 
finite densities the mode coupling approximation Eq. (|2^) can be used 

6^ 



a : 



lim Do 



$fc(w) • 



I). For 



(28) 



/o Sk 

In the limit of a macroscopic tracer, which becomes impenetrable to the Brownian 
particles so that /i'* (r) approaches a step function as argued in Eq. (^7|) , the integration 
in Eq. (28) already converges for fca <C 1. The density correlator is then given by the 
hydrodynamic limit and the structure factor equals the compressibility, Sk = Sq, so 
that 



(29) 



The leading contribution 2/3 arises from the static term, while we may again expect 
[^6| the Green-Kubo expression to sum up to the missing 1/3 if, extending our pair- 
density factorization, higher order density fluctuations are included. 



4. Application to diffusive long time tail 

An immediate consequence of the long-ranged structure built by particle diffusion 
around a macroscopic tracer, are slow time dependent fluctuations in the force 
the tracer feels. An interesting aspect of these so-called long time tails is that 
hydrodynamic calculations provide insig hts into them |,|o), even in the presence 
of boundaries |^l|] . As an application of the above discussion of boundary conditions 
for the diffusion equation, we study the long time tail in the force autocorrelation 
function of a tracer diffusing among Brownian particles. This extends the knowledge 



available at infinite dilution [p|, |18|,|19 



As a first step, the constitutive equation connecting the force a particle feels to 
the fluctuations of the conserved variable, the bath density, is required. It follows from 



the Zwanzig-Mori decomposition as |10 



(^F-(c.)) = {{FXe^} + (F-i?'(c^)f]^?q)) -> ^q-|(5e;(u;)) ,(30) 

where the limit in the second part holds for g — )■ and a; — )■ 0, and the coefficient a 
was defined in Eq. (^6|). Therefore, in the hydrodynamic limit, if a density gradient 
exists, it causes the force field: 

{SFs{r,u))^^V {Sg{r,co)) . (31) 

Do 

A sphere among the Brownian particles experiences this force density, and if 
it moves with velocity v, the density fluctuations in its vicinity are described by 
Eqs. (^,^6|). Following, the macroscopic approach to long time tails, and inserting 
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these expressions (with obvious definition of the matrix field a{r,uj)) into Eq. ( pT| ) 
gives the force field around the sphere. The sphere feels the interactions on its surface 
and thus the total force on it is obtained by averaging over the surface: 

F^,^ = -^v„ y -j^i"^^ ■ Vr) J d^svuj ■ a{s,uj) ■ $(r- s,uj) 



where $(r,w) = with — —iuj/D follows from Eq. ([T^). In general this 

result can not be simplified and e.g. the lu = value, which would be connected to 
the tracer diffusion coefficient, cannot be found from our hydrodynamic consideration 
alone. Linear response theory enables one to identify the tracer force autocorrelation 
function from Eq. (^: Fs^^j = — (l/3)( Fs('^) ' Fs) v^j. In the dilute hmit, it shows 
a small frequency anomaly of order and expecting this result at finite densities 
also, we expand the fundamental solution $ of the diffusion equation up to this order. 
Abbreviating the uninteresting terms this leads to: 

,,21 3/2 



47rL>S'ocr3 



D 



OK). (33) 



Importantly, the linear term in k vanishes (it would indicate {Fs{t — > oo) • Fg) cx 
t~^^^), and the leading w"^/^ anomaly corresponds to the final power law decay 
{Fs{t ^ oo) • F,) - 37rnaV[16S'o(7ri:)i)5/2]. As was expected from the spatial 
long-ranged pattern around the tracer there exists a temporal long time tail whose 
amplitude is closely connected to the latter. Interestingly, standard mode coupling 
theory gives different results ||,|l3|,|l^, ^ and requires additional considerations | p3[ |. 
Our result implies that a frequency dependent vertex coupling density modes would be 
required in order to obtain Eq. (^), which agrees with the known low density results 
p8| , upon accounting for the tracer diffusion by the replacement D ^ Dq + and 
the identification a = a+a^. The tracer mean squared displacement is connected to the 
force correlation function via the equation of motion, df{Ar^{t)) ~ —2Dl{Fg{t) ■ Fg), 
and exhibits a power law approach oc — a^t"^/^ to the long time diffusion. 



5. Conclusions and outlook 



We have presented the first statistical mechanics derivation of a hydrodynamic 
boundary condition for the diffusion equation, including the definition of the 
macroscopic boundary position from molecular parameters. This has proven 
surprisingly difficult, because a non-analytic (non-Markovian) structure in the 
resulting relaxation kernel needed to be identified. Quoting exact low density results 
and performing a mode coupling approximation, this structure could be established. 
The 1 / y/uj anomaly of diffusion close to a wall is connected to the penetration length 
y/D/to which arises generically in these situations Q. Thus, we believe the appearence 
of non-Markovian relaxation kernels is inherent to the statistical mechanics derivation 
of boundary conditions. Gratifyingly, within the mode coupling approximation it 
arises from a non-linear coupling of the hydrodynamic modes themselves, and thus, 
as required for a macroscopic hydrodynamic concept, does not depend on molecular 
details. 

As an application of the derived boundary condition, we calculated the long range 
density pattern around, the resulting force on, and the mean squared displacement of 
a tracer sphere moving among the bath particles. Extending results for dilute sytems, 
a generic power law approach (long time tail) to the ultimate diffusion was found in 
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the motion of the tracer, whose amphtude (cx a^) we determined for arbitrary tracer 
size. Interestingly, a may vanish for special interaction parameters, in which case the 
tracer might be called "invisible" . 

Clearly our calculation only presents a first step to determining boundary 
conditions for hydrodynamic equations. Thus when considering dispersed particles, 
their hydrodynamic interactions mediated via the solvent should be included. As these 
are long-ranged, interesting effects may appear. Also the description of the surface 
should be improved by going beyond the single tracer calculation. Structured or rough 
surfaces could be modeled, as well as fluctuating ones in order to adress the dynamics 
close to e.g. membranes. As the Smoluchowski equation is the basis for the dynamics 
of polymers and general macromolecules, the diffusion of complex molecules close to 
surfaces also could be analyzed following the presented approach. We hope to address 
some of these points in the future. 



Acknowledgments 

We want to thank M.E. Gates for his support and many helpful discussions and his 
critical reading of the manuscript. We also thank L. Bocquet and J.-L. Barrat for 
valuable discussions. M.F. was supported by the Deutsche Forschungsgemeinschaft 
under grant Fu 309/3 and through the SFB 563, and K.K. by a European Community 
Marie Curie Fellowship. 



Appendix A. Solution of the diffusion equation around a sphere 

For convenience of comparison with the results obtained in the main text, this abstract 
summarizes some hydrodynamic expressions. To linear order in the perturbing tracer 
velocity v, it is straightforward to solve the macroscopic Eqs. (|l|,H) around a sphere 
with radius a that oscillates with velocity amplitude v^^. With = —iui/D the 
square of the inverse skin depth, one obtains 

nu.{r) = • V ' ^ e--^--") . A.l 

2 + 2n<7 + [K,aY 

It simplifies to the far-field expression 

, , ^na^l + nr 2nina^ ^, 

n^{r) ^ v^-r -— —e , nq,„~v^-q— — , A.2) 

2D Dq^ — luo 

for a small sphere, cr ^ 0. Close to the sphere, it takes the near-field expression 



Mx)^-yI , n^^^^^uYl ^4 , A.3 

D K Dq^ — luj 

because the sphere degenerates to a plane upon taking the limit cr — s- cx); here x 
gives the distance ( r = cr x + r) and the Fourier transformation in Eq. ( [A.3| ) is one- 
dimensional. 



Appendix B. Mobile tracer calculation 

To simplify the presentation in the main text and to make direct contact with the 
hydrodynamic calculation, we have worked with a non-fluctuating, macroscopic tracer 
throughout. This appendix extends the calculations of the main text to a finite tracer 
mobility D^, and ckecks that they are indeed recovered by taking the limit Dg 0. 
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Concerning the tracer velocity in Eq. (^; with Eqs. (Q,^, and working in hnear 
response approximation, it becomes 

(9tr,)("<=) = (v) + (i^,F,)("^) = v-Dsf dr v(t).(F, e^^^*-^) F,) .(B.l) 

The deviation from the hydrodynamic velocity v is explicitly of the order Ds and thus 
vanishes as required for = 0. 

At finite D^, the resolvents i?q and R in Section 2.4 differ. Further, the tagged 
particle density now being a dynamic flucutating varible, it has to be considered as a 
separate slow mode in the projection in Eqs. ( pl] , p^ : 

P - Q%){q% + - ~ Q^)^{q% ■ (B.2) 

Here, Cq = ftq/S'q is the (tagged) direct correlation function, ft is straightforward 
to check that PP = P, Pgq) = Qq), and Pg^q) = f?q) up to corrections that are 
smaller by (n/iq)^/(7V5q) relative to the leading order. The susceptibility of Eq. ( p| ) 
is rewritten as 

nX^H = - [{F%Qq) - {F^R'^ilqgq)] {Q*^Rq0q)/{NSq) . (B.3) 

Hence, exactly the same decomposition of the susceptibility into an instantaneous and 
a retarded contribution has been achieved as in Eq. dl4) of the main text because of 
the negligible feedback of the tracer onto the bulk |jl^. To make contact with the 
expressions in the main text, a factorization approximation is required because of the 
tracer motion 

{e*qR<ig^) = {gqg%Rgqg^) ~ NSq<i>^{u;)<^>%{Lj) , (B.4) 

where $q(w) = {g"* R{lj) g"^) . The memory function ( F^q i?qf2qgq) also acquires a 
contribution from the tracer diffusion. To relevant lowest order in q it becomes 

*Do(F-i?'qq.Fq)-zD,(F:;i?'qq.F,eq) . (B.5) 

The limit — > recovers the results in the main text, where Dg — from the outset. 
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